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Abstract. In the paper [Fendley et al., J. Phys. A: Math. Gen., 38 (2005), 
pp. 315-322], Fendley, Schoutens and van Eerten studied the hard square 
model at negative activity. They found analytical and numerical evidence 
that the eigenvalues of the transfer matrix with periodic boundary were 
all roots of unity. They also conjectured that for an m x n square grid, 
with doubly periodic boundary, the partition function is equal to 1 when 
m and n are relatively prime. These conjectures were proven in [Jonsson, 
Electronic J. Combin., 13(1) (2006), R67]. There, it was also noted that 
the cylindrical case seemed to have interesting properties when the cir- 
cumference of the cylinder is odd. In particular, when 3 is a divisor of 
both the circumference and the width of the cylinder minus 1, the parti- 
tion function is -2. Otherwise, it is equal to 1. In this paper, we investigate 
the hard square and hard hexagon models at activity -1, with single peri- 
odic boundary, i.e, cylindrical identifications, using both topological and 
combinatorial techniques. We compute the homology groups of the associ- 
ated independence complex for small sizes and suggest a matching which, 
we believe, with further analysis could help solve the conjecture. We also 
briefly review a technique recently described by Bousquet-Melou, Linusson 
and Nevo, for determining some of the eigenvalues of the transfer matrix 
of the hard square model with cylindrical identification using a related, 
but more easily analysed model. 



1 Introduction 

Let G = {V, E) be an undirected graph. An independent set I in G is a subset of 
V(G) such that u, v € V(G) implies (u,v) £ E(G). Now, take G to be the graph 
in Figure [TJ where the half-edges are used to indicate that (1, 1) and (1,3) are 
adjacent and that (2, 1) and (2, 3) are adjacent. 

We are interested in the independent sets of G and start by listing them. We 
begin with the smallest set, the empty set 0, which of course is an independent 
set in any graph. Then, the singleton sets {(1,1)}, {(1,2)}, {(2,3)} are all 
independent. Finally, for each of the vertices (1, 1), (1, 2) and (1, 3) we can choose 
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Fig. 1. An undirected graph G. 



from two vertices of (2, 1), (2, 2) and (2, 3) to form a total of 6 independent sets 
of cardinality 2. We note that the alternating sum over all independent sets a, 

£(-i) kl = i- (i) 

cr 

It turns out that this happens a lot for graphs similar to G. Let the vertex set 
of C m ^ n be (i, j) for 1 < i < m and 1 < j < n and let there be an edge between 
and (i 2 , 32) if h = *2 and \j x - j 2 \ = 1 or if [*i - i 2 \ = 1 and jt = j 2 . 
Finally, let there be edges between each pair of vertices (i, 1) and (i,n). With this 
construction, we have G = 62,3. It is conjectured [5] that when m is odd and 3 
is not a divisor of gcd(m — 1,1%), then the sum ([TJ) over all independent sets a of 
C m ,n equals 1 and that if m is odd and 3 does divide gcd(m — 1, n), then this sum 
equals -2. 

A hard particle model in statistical mechanics is defined on an underlying 
graph. A configuration of the model is a set of particles placed on the vertices 
of this graph such that no two particles are adjacent. One imagines that the 
particles have a certain size, or shape, that extends to their neighbouring vertices 
in the graph. Furthermore, the particles may not overlap, as they are considered 
hard. On a square grid, one talks about the hard square model, because of the 
shape of the particles. On a hexagonal grid, the particles are triangles as each 
vertex have three neighbours. The correspondence between configurations of hard 
particle models and independent sets on the underlying graph is clear. 

The remarkable properties of the hard particle model on square grid graphs 
with cylindrical identifications were first observed by Fendley, Schoutens and van 
Eerten in 0]. They also studied the hard square model with doubly periodic 
boundary. Such identifications have the advantage that they make the graph 
vertex-transitive, i.e, the graph looks the same from every vertex. In this way 
the influence from the boundary is removed and these graphs are in that sense 
closer to the infinite square grid than the graphs obtained by simply taking an 
m by n subgraph of this grid. Fendley et al. conjectured that for the toroidal 
identifications, the roots of the characteristic polynomial of the associated trans- 
fer matrix were all roots of unity. Furthermore, they conjectured that when m 
and n were taken to be relatively prime, the alternating sum (fTJ) vanished. These 
conjectures were proven in 2006 by Jonsson [9] in a purely combinatorial way. 

To determine the partition function of these models at activity -1 has a direct 
correspondence in determining the Eulcr characteristic of the independence com- 
plex of the underlying graph. The matchings used in [5] were however not suitable 
for further topological analysis of the independence complex. Bousquet-Melou, Li- 
nusson and Nevo [3] studied a different region of the square grid which turned 
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out to be more easily handled using Morse matchings. They also determined the 
homotopy type of the independence complex on a parallelogram in the square 
grid and related these findings to certain eigenvalues of the transfer matrix for 
the hard square model on C TOi „. 

In this paper, we investigate the independence complexes for the square and 
hexagonal grids with cylindrical identifications. We obtain some results for small 
sizes by explicitly calculating the homology groups using a computer software 
package called Polymake [10] . We also construct Morse matchings for some classes 
of complexes with small circumference. The results obtained are enough to make 
certain conjectures on the homology groups. In particular, for odd n, there seems 
to be a symmetry in the square grid case and we conjecture that for i > and 
j, k > 1, we have 

ffi(/(Cj,2fc+l)) ^ Hi(I(C k , 2j +l)). 

However, it is also clear that the homology does not behave nicely in many cases, 
even for odd n, which clearly reduces the applicability of this approach. 

The paper is organised as follows. In Section [5] we present some basic notions 
and results from algebraic topology, discrete Morse theory and hard particle mod- 
els. We use a particular type of Morse matching to determine the homology of 
some independence complexes related to the hard square and hard triangle mod- 
els in Section [3] In Section [4] we introduce a different matching on the face poset 
of these independence complexes, which, while not being a Morse matching, still 
turns out to give some information on the alternating sum of the independence 
complexes of C m%n . Finally, in Section [5] we review an idea from [3] which exploits 
the knowledge of the partition function of a hard square model on a different 
region in order to try and obtain parts of the spectrum of the transfer matrix for 

2 Preliminaries 

In this section we define simplicial complexes, and in particular independence 
complexes on grid graphs. We also review the basic results from discrete Morse 
theory and give a framework for constructing Morse matchings on independence 
complexes. Finally, we talk about the connection between hard particle models 
from statistical mechanics and independence complexes on grid graphs. For a 
detailed introduction to algebraic topology, simplicial complexes and homology, 
see for example [T3] by Munkres. A classic reference for models in statistical 
mechanics is [2]. 

Throughout the paper we will use the following notation for addition and 
deletion of a single element x to or from a set S. 

S + x := SU{x}, S-x:=S\{x}. 

To be precise about our base cases, we define the Fibonacci numbers Fi in the 
following way. 

F = 0, Fi = 1, and F l+2 = F i+1 + F, for % > 0. 
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2.1 Simplicial Complexes 

Let X and Y be two topological spaces. If there is a continuous map / : X — > Y, 
with a continuous inverse, then X and Y are said to be homeomorphic, denoted 
by X = Y. 

Two continuous maps / : X — > Y and g : X — » Y are called homotopic if 
there is a continuous map F : X x [0, 1] — > Y such that 0) = /(x) and 
-F(;r, 1) = <?(a;). The spaces X and Y are said to be homotopy equivalent, denoted 
by X ~ Y, if there are maps / : X — > Y and /i : Y — > X such that / o h is 
homotopic to the identity map on Y and h o f is homotopic to the identity map 
on X. A space X that is homotopy equivalent to a point is called contractible. 
This will be denoted by X ~ •. 

An (abstract) simplicial complex A, is a collection of subsets over a ground set 
V such that if er <E A and r C a, then t £ A. The elements of Z\ are called faces, 
and the maximal faces are called facets. The dimension of a face tr is defined to 
be dim(er) = |er| — 1. 

The reduced Euler characteristic x of a simplicial complex A is the following 
sum. 

X(A):=J2(-lf imM (2) 

Given a simplicial complex A, the homology group Hi(A) of A in dimension 
i, over a ring R, is defined by Hi(A) — d(A)/ Bi(A). Here, the groups C, and 
Bi are the groups of i-chains and i-boundaries, respectively. The rank of Hi(A) 
is called the ith Betti number and is denoted by The following relation to the 
Euler characteristic, which can alternatively be taken as its definition, is called 
the Euler- Poincare formula. 

x(A) = J2 Mm (3) 

i>-l 

Here, /3_i is defined to be 1 if the complex A = {0} and otherwise. Homotopy 
equivalent spaces have isomorphic homology groups so x is a topological invariant. 

For a topological space X, we let dX denote the boundary of X and int X = 
X \ dX the interior of X. The n-dimcnsional simplex A n is the full simplicial 
complex on a n + 1-point ground set, i.e, A n = 2 V , with \V\ = n + 1. The n- 
dimcnsional ball B n is the subspace of the n-dimcnsional Euclidean space R™ 
given by B" := {v £ R™ | ||u|| < 1}. The ^.-dimensional sphere S n is similarly 
given by S n := {v £ « n+1 \ \\v\\ = 1}. That is, S n = dB n+1 . Apparently, we 
have B n = A n and S n = dA n+ i, where the simplicial complex dA n+ i can be 
described by 2 V \ V, for some |V| = n + 1. 

Let A be a simplicial complex. The cone over A, denoted by conc(Z\), is the 
simplicial complex over the ground set of A together with a new vertex c, the 
cone point. 

cone(Z\) = {er, <t + c | a £ A} 

The cone over any simplicial complex is contractible. 

The suspension over A, denoted by susp(Z\), is the following simplicial complex 
over the ground set of A together with two new vertices sq and s\. 

susp(Z\) = {a, a + s , a + s\ \ a £ A} 
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For any simplicial complex A and i > 0, 



Hi(A) S H i+X (susp(Zi)) . (4) 
2.2 Independence Complexes of Regular Grids 

Regular grids arc tilings of the 2-dimcnsional plane by regular convex polygons. 
There exist three types of regular tilings, triangular, square and hexagonal. In this 
paper we have considered the square and the hexagonal grid. These are depicted 
in Figure [2 




Fig. 2. (a) The square grid, (b) The hexagonal grid. 



Let S be the graph induced by the boundary of the square grid. Then, S 
is isomorphic to an infinite graph on the vertex set Z x Z, with edges between 
(til, U2) and (vi, V2) if u\ = v\ and \ui — v<i\ = 1 or if \ui — Vi\ = 1 and u-i = V2- 
From now on, we will assume that S is defined in this way. Let <SVn,n be the graph 
induced from S by the vertex set [m] x [n], where [m] denotes the set {1, . . . , m}. 
By identifying sides of S m ^ n , we can obtain graphs embeddable on different 2- 
dimcnsional topological surfaces. Specifically, let G m . n be S m , n +\, with vertices 
(i,n + 1) and (£, 1) identified, for i G [m]. Furthermore, let T m%n as the graph 
'S'm+ijn+ii with the vertices (i, n + 1) and (i, 1) as well as the vertices (m+1, j) and 
(1, j) identified, for i G [m] and j G [n], respectively. The graph C mj „ is obviously 
embeddable on a cylinder (but not on the plane), while T mj „ is embeddable on a 
torus (but not on a cylinder). 

Let H be the graph induced by the boundary of the hexagonal grid. Similarly 
to the case of S 1 , we can find an appropriate representation of the graph H with 
vertex set Z x Z. This graph has edges between (ui,U2) and (ui,i>2) if ui = vi 
and \u2 — ua| = 1 or if u\ — v\ — 1, U2 = V2, and u\ + =2 0. For positive 
integers m and n let H m .n be the subgraph of H shown in Figure E^a). We can 
straighten this graph and draw it on the square grid, as is shown in (b). We will 
use this brick wall-like appearance in this paper. If we identify the left and right 
boundaries of the graph in Figure [3f 6) , we get a graph n with the topology of 
a cylinder. Figure [3fc) shows an alternative way of drawing this cylinder, using 
half edges on the left and right boundaries to indicate the identifications. We will 
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(7,8) 



(a) 



(1,1) (1,8) 
(c) 



Fig. 3. (a) H 6: 4 as a subgraph of H. (b) iJ 6j 4 as a brick wall, (c) The cylinder 
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use the vertex set [m] x [n] C Z x Z to represent these graphs, with the vertices 
numbered as in Figure [3tc) ■ We will also briefly look at the hexagonal grid with 
toroidal identifications. In this case, we will use the representation in Figure^fo), 
identifying the left and right boundaries as well as the top and bottom ones. We 
denote the graph obtained by these identifications by „. Note that this graph is 
different from the one we would obtain by doing the top to bottom identifications 
in Figure [3Jc). 

We will be studying simplicial complexes related to these regular grids. Let 
G = (V, E) be a finite, simple and undirected graph. An independent set I in G 
is a subset of the vertex set V(G) such that if u, v € I then (u, v) $ E(G). That 
is, no pair of neighbours may be in / at the same time. 

Define the independence complex 1(G) as the simplicial complex on the ground 
set V(G) with a g 1(G) if and only if a is an independent set in G. Similarly, 
let X(G) be the flag complex of G. This is the simplicial complex on the ground 
set V(G) with a 6 X(G) if and only if a induces a complete subgraph in G. It 
is obvious from the definitions that 1(G) = X(G), where G = (V, E) denotes the 
complement graph of G with respect to the complete graph on V, i.c, E = (]^)\E. 

In general, computing 1(G) is believed to be hard, since the facets of 1(G) 
are the maximal independent sets of G. Any such facet of maximal dimension 
therefore gives a solution to the maximum independent set problem, which is a 
well known NP-completc problem Dual to this is the NP-complete maximum 
clique problem of finding a complete subgraph (or clique) in G of maximal size. 
These cliques are given by the highest dimensional facets of X(G). 



2.3 Discrete Morse Theory 

Discrete Morse theory was introduced in [5] by Robin Forman. It relates the 
homotopy type of a simplicial complex to that of a CW-complex with (ideally) 
simpler structure. See [6] for a hands-on introduction to this subject. 

Let A be a simplicial complex. We define the face poset P(A) on A as the set 
of faces in A ordered by inclusion. The underlying graph of the Hasse diagram 
of P(A), or simply the underlying graph of P(A) 1 is the directed graph with the 
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faces of A as vertices and the cover relations of P(A) as edges, directed from the 
smaller faces to the larger faces. Let M be a matching on the underlying graph of 
P{A). We will also say that M is a matching on the complex A. The graph induced 
by M is the graph obtained from the underlying graph of the Hasse diagram by 
reversing the direction of the edges which are in M. 

Definition 2.1. A Morse matching M on a simplicial complex A is a (partial) 
matching of the Hasse diagram of P(A) such that the graph induced by M is 
acyclic. 

The main theorem of Morse theory relates the homotopy type of a simplicial 
complex to that of a certain CW-complex related to the unmatched faces of a 
Morse matching. These unmatched faces are sometimes referred to as the critical 
cells of the matching. For our purposes, the following definition of a finite CW- 
complex will suffice. 

Definition 2.2. A finite CW-complex (or cell complex) C is a collection of cells 
such that for each cell e a 

(a) There is a homeomorphism from B k to e a for some k. 

(b) e a \ int e a is the union of cells eg in C . 

The main theorem of discrete Morse theory can be stated as follows, using the 
notion of a Morse matching. 

Theorem 2.1. Let A be a simplicial complex with a Morse matching M . Assume 
that for each p > 0, there are u p unmatched p-simplices. Then A is homotopy 
equivalent to a CW-complex with exactly u p cells of each dimension p, for each 
p > 0, and uq + 1 cells of dimension 0. 

Corollary 2.2 Let A be a simplicial complex with a Morse matching M such 
that Up = for all but one p. Then, for this particular p, 

A~\/S p . 

Up 

That is, A is homotopy equivalent to a wedge of u p p-dimensional spheres. 

Corollary 2.3 Let A be a simplicial complex with a perfect Morse matching M . 
Then, A is contractible. 

In this paper, we will apply discrete Morse theory to independence complexes. 
The acyclic matchings that we will use are all instances of the following general 
construction. By describing them in this framework, we can easily determine their 
acyclicity. The construction is similar to the matching trees of [3] . Throughout this 
paper, we will however use the term matching tree to refer to our construction. 

Example 2.1. The following example shows the basic idea behind the matching 
trees. In Figure |4j we construct a matching on the complex /(S3.2). First we 
partition the set of faces in to two disjoint sets, based on whether the vertex (1, 1) 
is in the face or not. We indicate this by writing (1,1) to the left of the root in 
the tree. To the right of the root we place the graph S^,2 which indicates that 
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all faces of the complex are still unmatched. In the left subtree we have forced 
(1, 1) not to appear in any face. We illustrate this by placing an empty particle 
on (1, 1) in the graph to the right of the root of this subtree. In the right subtree 
we have instead forced (1,1) to be present in every face. This is illustrated by 
a filled particle on (1,1) and empty particles on its neighbours (1,2) and (2,1). 
In both subtrees, this effectively reduces the problem to the construction of a 
matching on two smaller graphs, induced from 5^2 by the complement of the 
fixed vertices. In the next step, the vertex (1,2) in the left graph, and (2,2) in 
the right graph both have degree 1. We use these vertices to create the matchings 
M(l,2) and M(2,2) in the respective subtree. Let E t = I(S 3 . 2 - (1,1)) and 
S r = I(Ss,2 \ {(1, 1), (1, 2), (2, 1)}), and define the matchings as follows. 

M(l,2) :={(*,*- (1,2)) \ (1,2)6^6^}, 

M (2, 2) := {(a, a - (2, 2)) | (2, 2) G a G E r }. 

The unmatched faces a remaining after removing the faces matched by M(l, 2) 
must have (1, 2), (2, 1), (3, 2) g" a and (2, 2) G a. This further reduces the graph. In 
the right subtree, it completely fixes one single graph. In the left subtree, there is 
still a vertex (3, 1) on which we can match. Matching {(2, 2)} with {(2, 2), (3, 1)} 
removes all the faces from this subtree, which we indicate by writing to the right 
of the last node. 



(3,1) (3,2) 




Fig. 4. A matching tree for the simplicial complex liS-sa)- 



In Proposition I2.6[ we wm see that the union of the matchings made in this 
type of construction forms an acyclic matching on the original complex. We can 
therefore apply Theorem 12.11 and due to the remaining 1-dimensional face in the 
right subtree, which is the only critical cell, we conclude that 1(53.2) is homotopy 
equivalent to a 1-dimensional sphere. 
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Definition 2.3 (Matching tree). Let T be a finite, rooted tree in which each 
internal node has either one or two children. To each internal node (which includes 
the root) t G T, we associate a vertex p(t) G G. 

We say that T is a matching tree if the set of vertices p(t) fulfil the following 
two conditions. For any node t r G T, let (ti, . . . ,t r ) be the path from the root t\ 
to t r . Then, 

— all p(ti), . . . ,p(t r ) are distinct, and 

— {p(t) | t G {ii, . . . , t r } H V m } forms an independent set in G, 

where V m denotes the set of nodes with one child. 

The vertices p(t) are called the pivots in [3]. They are the vertices that we 
wrote to the left of the nodes in Figure 0J Given the complex 1(G) and the tree 
T, we determine a matching on P(I(G)) as follows. To each node t G T, we 
associate a content C(t), which is a subset of the complex 1(G). This is the set 
that was illustrated by the graphs to the right of the nodes in Figured] For the 
root r G T, we let C(r) = 1(G), and for the other nodes the content will be 
determined recursively. 

If t has two children, then we call t a split node, and we let the content of 
the root of the left subtree of t be the set of faces {cr G C(t) \ p(t) g" a} and the 
content of the root of the right subtree of t be the faces {a eC(i) | p(t) G cr}. 

If t has only one child, then we call t a match node. We define the partial 
matching M(t) associated to the match node and the set of matched faces F(t) 
as follows. 

M(t):={(a,a~p(t))\p(t)eaeC(t)}, F(t) := |J {a,r} 

(a,T)eM(t) 

We will see below that the set M(t) forms a partial matching on C(t). For match 
nodes, we let the content of the child of t be C(t) \ F(t). 

Definition 2.4. Let G be a graph and a G 1(G). We say that a position x G V(G) 
is free in a if x G" a and a + x is an independent set in G. 

The following lemma provides all the structure we need on the sets C(t) to show 
acyclicity for our matchings. 

Lemma 2.4 Let T be a matching tree and let t G V m . If a G C(i) and p(t) 
is free in a, then a + p(t) G C(t) Conversely, if a G C(t) and p(t) G a, then 
<J-p(t) G C(t). 

Proof. Both parts of the lemma obviously hold true for the root node. We will 
show that if either part fails for some t, where t is taken to be the first node on 
the path from the root to t for which this is the case, then the other part must 
fail at an earlier node t' , closer to the root. If we apply this argument twice, we 
reach a node t" which is strictly between the root and t and for which the same 
part of the lemma fails. This contradicts our choice of t, and the lemma follows. 

Assume that a G C(t) and that p(t) is free in a. The reason for a + p(t) not 
to be in C(t) could be that there was a split node which had a in one subtree 
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and <r+p(t) in the other. This can not happen, however, since the two faces only 
differ in p(t) and along the path from the root, all pivots must be distinct. 

The only remaining possibility is then that a +p{t) was matched in an earlier 
node t' . Now, if p(t') is free in a + pit), then it is also free in a, and by our 
assumption, a would have been matched at t' as well. So, we can assume that 
pit') £ cr + pit), and then also pit') £ a. The second part of the lemma therefore 
fails at t'. 

For the second part, assume that a £ Cit) and p(t) £ a, but that a — pit) £ 
C(t). Again, a — pit) could not have been removed by an earlier split node since 
it only differs in p(t) from a. 

So assume that a — pit) was matched in an earlier node t' . If pit') £ a — pit), 
then pit') £ a and by induction, a would have been matched with a — pit') at t' . 
Therefore, pit!) is free in a — pit). By the restriction that the pivots of the match 
nodes form an independent set in G, we have that pit!) is free in a as well. The 
first part of the lemma therefore fails at t'. □ 

Using Lemma E3J it is easy to see that in fact, for t £ V m , we have 

Fit) C Cit). 

That is, for each match node t £ V m , M(t) is a partial matching on the subgraph 
of the underlying graph of _P(/(G)) induced by the set Cit). It is acyclic, as the 
only way to "move down" in the directed graph induced by M(t) is by removing 
pit). 

It therefore makes sense to define the matching M determined by the tree T 

as 

M := |J M(t). (5) 

tev m 

We will now show that M is acyclic, so that Thcorem l2.1l is applicable. The critical 
cells will be those in the set Utev; where Vi denotes the set of leaves in T. 
Our proof will use the following lemma, which appears in [S] . 

Lemma 2.5 (Cluster Lemma) Let A be a simplicial complex and f : Pi A) — > 

Q be a poset map, where Q is an arbitrary poset. For q £ Q, let M q be an acyclic 
matching on f~ 1 iq). Let 

M= (J M q . 

Then, M is an acyclic matching on A. 

The following proposition justifies our definition of a matching tree, show- 
ing that any tree gives an acyclic matching on the corresponding independence 
complex. The proof is entirely guided by the Cluster lemma (Lemma I2.5|) . 

Proposition 2.6 Let G be a graph and T a matching tree defined on the inde- 
pendence complex /(G). Then, the matching M given by (fj| on /(G) is acyclic. 

Proof. We start by defining a total order on the nodes of the tree T. First, if t' is 
an ancestor of t, then let t <t t'. Then, if T\ and T2 are two subtrees of a node 
t, with T\ being the subtree where pit) £" a for any a £ Crti), let t\ <t ti for all 
t x £ T x and t 2 £ T 2 . 
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Let Q = (V(T),< T ) and 



r\t) 



F(t) ifteV m , 
C(t) iSteVu 

otherwise. 



It is easy to check that the F(t) are pairwise disjoint and that F(ti) Cfa) = 
for all ti G V m and G Vj. Furthermore, each face of /(G) is either in F(t±) for 
some t\ G K« or in C(<2) for some £2 G VJ. Thus, / : /(G) — * Q is a uniquely 
defined map. It remains to verify that it is indeed a poset map. 

Let <7, t G /(G), with a C r. If /(cr) and /(r) are in different subtrees of some 
node f , then a must be in the left, as we know that p{t) is in all faces of the right 
subtree, but in none of the left. Thus, f(a) <t /(t). 

If instead /(cr) and /(r) are on the same path from the root to some leaf, then 
we will assume that /(r) is closer to the root than f(cr). Let v = p(f(r)) be the 
pivot at /(t). If v G cr, then Lemma [2~4l says that cr should have been matched in 
the node /(t), so we have a contradiction. Therefore, u ^ cr. But then a C r — w, 
and t — v was matched with r. This implies that v is free in cr, and Lemma 12.41 
again shows that a should have been matched in /(t). We have a contradiction. 
Therefore f(r) must be equal to or lie further from the root than /(cr), which 
means that /(cr) < T /(r). 

Finally, for each t G V m , we have that M(t) is a perfect, acyclic matching of 
/ _1 (i) = F(t). We can therefore use Lemma [231 to conclude that M is acyclic on 
the entire /(G). □ 

We finally note that our construction disallows some strategies for choosing 
split and match nodes which would have resulted in acyclic matchings. This is 
due to the constraint that the match nodes in each path from the root must form 
an independent set. This condition is not strictly necessary, but for the purposes 
of this article, this construction will suffice. 

2.4 Hard particle models 

A particular type of model in two-dimensional statistical mechanics is the hard 
particle model. In such a model, given a regular grid like the ones described in 
Section 12.21 a configuration is a set of particles placed on the sites of the grid 
so that no two particles are adjacent. The name is explained by replacing each 
particle with a polygon, using the neighbours of the original particle as the vertices 
of the polygon. This polygon is assumed to be "hard" , meaning that no two 
polygons may overlap in a configuration. When the grid is the square grid, the 
polygons are squares, and the model referred to as the hard square model. For the 
hexagonal grid, the polygons are triangles, and one has the hard triangle model. 
For a thorough introduction to models in statistical mechanics, see [2], in which 
can also be found the solution to the hard particle model on a triangular grid. 
This solution was first presented by Baxter in [1] in 1980. 

To a model, one can associate a partition function Z. This is the generat- 
ing function of the configurations with respect to the number of particles of a 
configuration. 




(6) 



a 
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where the sum is taken over all allowed configurations of the model. The variable 
z of Z is called the activity and controls the behaviour of the model. 

A commonly used framework for deriving the partition function is the transfer 
matrix method. See [TH Chapter 4.7] for an introduction to this method. Here, 
we illustrate the idea by giving a transfer matrix for the rectangle S m ^ n . More 
examples are given in Section [5] 

Let L n denote the set of independent sets on the line Si, n . Since \Li\ = 2, 
|L 2 | = 3, and the recursive relation |£ n +2| = |£n+i| + \L n \ holds, we have \L n \ = 
F n+ 2, the Fibonacci numbers. Let T(z) := T(z) aT be the F n+2 x F n+2 matrix 
where the individual matrix elements are given by 



T{z) c 



Z \ T \ ifo-nT = t 

otherwise. 



The matrix T(z) is called a transfer matrix of the model. It describes the effect 
of adding an additional row. We will only be studying models at activity -1 so 
we leave out the variable z and write T = T{— 1). The generating function of 
Z(S mjl ) with respect to n can now be written as 

J2Z(S m , n )t n = (l-tT)- 1 . 
n>0 

In addition, by simply taking the trace of the nth power of the matrix T, we can 
also express the partition function Z(C m ^ n ) of a cylinder. 

Z(C m , n ) = tr(T(z)"), 

or as a generating function 

J2 Z (Cm,n)t n =tr(l -IT)- 1 . 

n>0 

The characteristic polynomial P(t) of T is given by 

P(t) = det(iJ - T). 

The roots of the characteristic polynomial can also be found by taking the poles 
of the generating function tr(l — tT) . 

It was observed by Fendley, Schoutens and van Eerten [1] that the hard square 
model (on a torus) exhibit some interesting properties when the activity is chosen 
to be —1. Among other things, they conjectured that all the roots of the charac- 
teristic polynomial arc roots of unity. These conjectures were rigorously proven by 
Jonsson [3] by means of applying an ingeniously designed series of matchings on 
the set of configurations. Doing this, Jonsson discovered an underlying structure 
of a particular type of rhombus tilings of the plane which he used to settle the 
conjectures. An attempt to mimic the idea of Jonsson is carried out in Section 21 
albeit with limited results. 

The following conjecture for the hard particle model on a cylinder in the 
square grid at activity -1 can be found in [3]. The conjecture is based on results 
for m < 11. 
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Conjecture 2.1. For odd n, 




1 



2 



if 3 | gcd(m — 1, n), 
otherwise. 



From @ and ^ we see that, at z 



-1, 



Z(G) 



X(I(G)). 



That is, the partition function at activity —1 is (up to a negative sign) equal to 
the reduced Euler characteristic of the independence complex of the underlying 
graph. For this reason, we will sometimes refer to Z(G) as the alternating sum 
of the corresponding independence complex. In the next section, we investigate if 
this connection to topology can be taken further than just to the reduced Euler 
characteristic. 

3 Explicit results 

In this section we present the explicitly calculated homology groups of some in- 
dependence complexes on graphs related to the square and the hexagonal grid. 
First, we do some small cases on the square grid by hand, for a fixed circumfer- 
ence n. Then, we present some larger cases, where we have used a computer to 
calculate the homology groups. The calculations were made using the software 
package Polymakc [TU] version 2.2 by Michael Joswig and Ewgenij Gawrilow with 
many contributors. 

3.1 Results for small circumference 

For some cases with small n, we can find suitable matching trees and work out 
the homotopy type. We do this here for the square grid when n < 5 and for the 
hexagonal grid when n = 1 or n = 2. For the case n = 4, we are only able to 
determine the alternating sum. This result is given in Proposition 13.71 We note 
that for the square grid, these results are consistent with Conjecture 12. II 

When n = 1 there can be no particles on the grid since every vertex in the 
grid is a neighbour to itself. Therefore /(C mj i) = {0} and Z(C m ,i) = 1. 

Proposition 3.1 Letm > 1. Then, I(C m>2 ) ^ S^ 2 ^ 1 andZ(C m - 2 ) = (-1)^/21-1. 

Proof. Note that C m ,2 = <Sm.2- The construction in Example 12.11 in Section [2] 
is easily extended to arbitrary m. The only unmatched configuration has to/2 
particles when m is even and (to + l)/2 particles when to is odd. □ 

Proposition 3.2 Let m > 1. Then, 



' 52m/3-i if to = 3 0, 

I(C m , 3 )~ \ V 2 S 2(m - 1)/3 ifm= 3 l, 
£2( m -2)/3+l ifm= 3 2, 



1.3 



(1,3) 



(2,1) 




Fig. 5. Matching tree for C t 



m , 3 



and 



1 ifm= 3 0, 
Z{C m ,a) = { -2 ifm = 3 1, 
1 ifm = 3 2. 



Proof. Let C' n n = C mtn — (1,1)- Wc construct an acyclic matching on C m . 3 as 
follows. If we start by splitting on (1,1), we are left with C' m 3 in the left subtree, 
and after a matching, a graph isomorphic to C' m _ 1 3 in the right subtree. Figure [3~Tl 
shows the matching tree construction for C' m 3 . Note that the last remaining graph 
in the left subtree is isomorphic to C' m _ 3 3 and that we have added two particles 
on the path from the root to this node. 

When m = 1, we can construct a matching leaving two configurations with 
one particle each, when m = 2, we find a perfect matching, and when m — 3, we 
are left with one configuration, having two particles. The proposition follows by 
adding the configurations from the two subtrees, and using Proposition 12.61 and 
Corollary O □ 



if m is even, 
if m is odd, 



Proposition 3.3 Let m > 1. Then, 




and 



m,4 ) 



m + 1 if m is even, 
—m ifm is odd. 
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Fig. 6. Matching tree for C mA - (1, 1). 



Proof. Let C' m n = C m>n — (1, 1). Figure l3~Tl shows how to construct a matching 
tree for the graph C' m 4 . Note that the left subtree gives a recursive relation of 
length 2. When m is even, there will be one configuration with m particles left 
at the end of this subtree. When m is odd, we have a perfect matching. In the 
right subtree, there remains a graph which is isomorphic to C' m _ 1 n . When m = 1, 
this subtree has one configuration with one particle. Thus, the entire tree leaves 
[m/2j + 1 configurations with m particles. 

Now, we build the matching tree for C m ^ as follows. We first split on (1, 1). 
The left subtree will be the one previously constructed for C' m4 with [m/2\ + 1 
remaining configurations. In the right subtree, we can match on (3, 1), leaving us 
with a graph isomorphic to C' m _ 2 4, and \_(m — 2)/2\ +1 remaining configurations. 

We then use Proposition 12.61 to show that we have an acyclic matching with 
2[m/2j + 1 unmatched configurations, each with m particles. Corollary 1 2 . 21 finally 
gives us the homotopy type. □ 

Proposition 3.4 Let m > 1. Then, 

if m is even, 
if m is odd, 

and 

Z(C m , 5 ) = 1. 

Proof. We give a sketch of the matching tree for C m ^. First, we split on both 
u = (1, 1) and v = (3, 1). For S C 2^^>, let 

X m (S) ={ae I(C m>5 ) I a n {u, v} = S}. 
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Then, X m ({u}) and X m ({v}) can be perfectly matched. If m is even, then after 
matching X m (0), one configuration with m particles remain, and there is a perfect 
matching for X m (%) otherwise. For X m ({u, v}), we have one configuration with 
77i + l particles left if m is odd and a perfect matching otherwise. □ 

Proposition 3.5 Let m > 1. Then, 

{ S 2m/3+i if m = 3 Q 

• if m = 3 1 

5 2(m+l)/3 if m = 3 2 



Fig. 7. The graph C4 1 is isomorphic to Pg . 



Proof. It can be seen in Figure [7] that the graph C^f j is isomorphic to P^i+i, the 
path on 2n + 2 vertices. The result immediately follows. □ 

Proposition 3.6 Let m > 1. Then, the simplicial complex I(C^ 2 ) * s homotopy 
equivalent to a wedge of F m +2 m- dimensional spheres. 

Proof. The Morse matching for I(C^ 2 ) ^ s given by the matching tree in Figure [5] 
We start by splitting on vertex v 2 . If there is a particle on v 2 , then we can find 
a perfect matching on all remaining configurations by matching on V4. If v 2 is 
empty we split on vertex U4 . The remaining graph in the right subtree is shown 
in the lower right corner of the figure. We see that it is isomorphic to C^_ t 2 , 
and that we have added one particle on the way to this node. In the left subtree, 
we eventually reach a node where the remaining graph is the one in the lower left 
corner. This graph is isomorphic to C^_ 2 2 and the two match nodes on the way 
to this node adds two particles. 

The base cases when m = 1 and m = 2 are constructed by hand. We can 
construct a matching on I(C^ 2 ) so that two configurations remain, both with 
two particles. On I{C 22 ) we can similarly construct a matching which leaves 
three configurations, each with three particles. Using these matchings and Propo- 
sition [223 we finally construct an acyclic matching on 7(C^ 2 ) which leaves F m+2 
unmatched configurations each with m + 1 particles. □ 

As we will see in the next section when we determine some homology groups 
using computer calculations, the complex L(C^ 4 ) has non-zero homology in two 
dimensions. Therefore, we can not in general hope to apply Corollary 12.21 to the 
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V3 




Fig. 8. A recursion for I(C^ 2 ). 



complexes /(C„ 4 ). However, Table \E\ suggests that the alternating sum may be 
simple. We are able to find a matching tree which shows that this alternating 
sum is identical to that of I{C^ 2 ) which was determined in Proposition 13. 61 This 
results is given in the following proposition. 

Proposition 3.7 Let m > 1. Then, Z(C% >4 ) = (-l) m+1 F m+2 . 

Proof. We give a sketch of the matching tree construction used to prove Proposi- 
tion !3.71 We start by splitting on the vertices V\, . . . , v± in Figured] The first thing 
to note is that if particles are present on two consecutive Vi and (modulo 
4), then we can match everything in this subtree. The upper right configuration 
shows this case, when for all faces a in this subtree, Ui,i>2 G &■ Here, we match 
on 115. A similar situation occurs when only one Vi is present. 

Therefore, we have either the case vi,Va € a,V2,V4 o (or the complement), 
or no Vi is in a. The lower left configuration shows the first case, where we proceed 
by matching on and vq. In this case, which occurs twice, we end up with the 
graph C^_ 2 4 and 4 particles fixed in each configuration. 

The lower right configuration shows the case when none of the Vi are in a. 
We then match on u±, . . . ,1*4. In the second case, we reach the graph C^_ 3j4 
after fixing 7 particles. We see here that the critical cells will be spread out over 
several dimensions. We therefore focus on the alternating sum, disregarding the 
homotopy type in this case. From the previous argument, we have the recursion 

Z{Cm,i) = 2Z(C%_ 2i i) - Z(C m _ 3 4 ), 

and base cases are 

Z{C* A ) = -1, Z(C? A ) = 2, Z{C? A ) = -3, 

where Cq 4 is isomorphic to an 8-cycle. The proposition follows by solving this 
recursion. □ 
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V5 «6 

Fig. 9. Split and match nodes in a matching tree for I(C m ^). 



3.2 The matching M 

In the rest of this section we will present the homology groups of some inde- 
pendence complexes which have been computer generated. To reduce the size of 
the complexes and the time of the calculations, we first apply a particular Morse 
matching to each complex. We start by introducing this matching, and the ma- 
nipulations made on the complex, before giving our results in Section 13.31 and 
Section l3~4l 

Let G be a simple, undirected graph. Let O € 1(G) be an independent set in 
G, which we will call the odd vertices. We will assume a fixed, but arbitrary total 
order < on O. 

The matching Mo on the Hasse diagram of P(I(G)) is defined as follows. 
The edge (a, a') is in Mo if o' \ a = {o}, where o G O is the smallest element 
(with respect to <) in O for which a U {o} is an independent set. According to 
Proposition 12.61 the matching Mo is acyclic, and therefore a Morse matching. 
According to Theorem 12.11 we then know that 1(G) is homotopy equivalent to 
some cell-complex, with cells in the dimensions given by the sizes of the unmatched 
faces in 1(G). We call this set of unmatched faces X and note that we can express 
it as 

X = {a G 1(G) | Vu G O : of] Nbd(u) ^ 0}, (7) 

where Nbd(u) denotes the neighbourhood of u in G. Of course, in general, X 
does not form a simplicial complex. From we see that the total order < does 
not affect X. Now, given a set O, and thereby a set X, we define two simplicial 
complexes Aq and Fo. 

Ao = {<T\<TCa' eX}, r = A \x 

We will now prove a theorem that relates the homotopy type of the original 
independence complex to that of Fq. First, however, we will give a topological 
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lemma. In [TO] it was stated for the case when A is a simplex. This slightly more 
general version is needed for the proof of our theorem. 

Lemma 3.8 Let r C A be simplicial complexes and assume that A is con- 
tractible. Then, A/T ~ susp(-T). 

Proof. The concept of a mapping cone is described in [7J Chapter 0]. In our 
case, the mapping cone is the simplicial complex X — A U;d conc(_T), i.e, we 
raise a cone in A over the subcomplex T. Note that X/cone(r) ~ A/T and 
that X/A ~ susp(_T). Now, since both cone(-T) and Z\ are contractiblc, we have 
X/cone(-T) ~ X ~ X/Zi which finishes the proof. □ 

We are now ready to state and prove the following theorem, which allows us 
to simplify our calculations of the homology groups. 

Theorem 3.9. Let O G 1(G) be an independent set and assume that Ao is 
contractible. Then, 

1(G) ~ susp(rb)- (8) 

Proof. The set X is given by J7J). Another way to express this is to say that a G X 
if and only if O C U ug(T Nbd(u), From the latter definition we can immediately 
see that if a C a' and a E X, then we have a' G X. Let Zo = 1(G) \ X. The 
previous argument shows that So is a simplicial complex. Figure [TO] illustrates 
the relationship between X and the simplicial complexes Aq , i~b an d ■ 






Fig. 10. The set X and the simplicial complexes Aq, 2~b an d 



Since Mo gives a perfect matching on w e can collapse So to a point in 
1(G). But this is the same as identifying 2~o with a point in zio- Thus, 

1(G) ~ J(G)/^o - ^o/^o- 
Using Lemma 13.81 we now see that /(G) ~ susp(Io). □ 
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The reason for using Jo instead of 1(G) is two-fold. First, we have used a 
matching to reduce the size of the complex, and thereby, hopefully, reduced the 
time to compute its homology. Secondly, in many cases a description of the com- 
plex Jo is easier to obtain and more compact than a description of the original 
complex. 

3.3 Graphs on the square grid 

Let G = C m , n and assume that n is even. In this case we choose the independent 
set O — {u = G V(G) | i+j is odd}. Then, Aq — an d thus 

contractible. The complex 2~o can be generated from the set O in the following 
way. For each u G O we have a maximal face {v G V(G) \ \ v £ Nbd(u)} G 
Jo- We see that for the maximal faces a G Jo, either dim(cr) = mn/2 — 5 or 
dim(cr) = mn/2 — 4, depending on whether the corresponding v G V(G) \ O is on 
the (non-identified) boundary of the cylinder or not. 

Example 3.1. Let m = 3, n = 4, and O = {(1, 2), (1, 4), (2, 1), (2, 3), (3, 2), (3, 4)}. 
Then, the six maximal faces of Jo are as shown in Figure [TT] and the complex 
itself can be realised as in Figure [T2J We can see that it is homotopy equivalent to 
a wedge of three 1-dimcnsional spheres. Therefore, we conclude from Thcorcm l3.9l 
that I(C 3A ) ~ V 3 S 2 . 



(3,1) (3,4) 




(1,1) (1,4) 



Fig. 11. The six maximal faces of Jo- 



For odd n, we let O = {u = G V(G) \ j ^ n,i +j is odd}. The complex 
Ao is still contractible but the maximal faces of Jo are less regular. 

Table [T] and Table [5] show the calculated homology groups for the complex 
I(C m ^ n ) for small m and n. An entry (k, d) in the table should be read as 
Hk(I(C m ^ n )) = Z d . Missing entries are cases where Polymake ran for too long, 
or more often, ran out of memory. The homology vanishes in all dimensions that 
are not listed. Note that for n < 5, the homology groups in these tables can be 
derived, for arbitrary m, from the results in Section 13.11 For odd n, we observe 
an interesting symmetry, which we pose as a conjecture. 

Conjecture 3.1. For all i > and j, k > 1, the following holds. 

Hi(I(Cj : 2k+i)) — Hi(I(Ck,2j+ij) 
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Fig. 12. The simplicial complex Jo — V3 S 1 - 



Table[3]is taken from [5] and shows the values of Z(C m ^ n ) for small m and n, 
calculated using transfer matrices. The cases when n is odd are left out as they 
all satisfy Conjecture 12. 11 In these cases, it was also observed, and conjectured to 
hold in general, that all the roots of the characteristic polynomials were roots of 
unity. In particular, this implies that the rows of Table [3] are periodic. 





m 




1 


2 


3 


4 


5 


6 


7 


8 




2 


(0,1) 


(0,1) 


(1,1) 


(1,1) 


(2,1) 


(2,1) 


(3,1) 


(3,1) 




4 


(0,1) 


(1,3) 


(2,3) 


(3,5) 


(4,5) 


(5,7) 


(6,7) 


(7,9) 




6 


(1,2) 


(2,1) 


(3,1) 


(5,4) 


(6,1) 


(7,1), (8,2) 


(9,4) 


(11,7) 


n 


8 


(2,1) 


(3,3) 


(5,5) 


(7,5) 


(8,1), (9,4) 


(11,7) 








10 


(2,1) 


(4,1) 


(7,1) 


(8,1), (9,2) 












12 


(3,2) 


(5,3) 


(8,3) 


(11,8) 












14 


(4,1) 


(6,1) 


(9,1) 













Table 1. Homology of I(C m . n ) for even n. 



Both for even and odd n, we find cases where there is non-zero homology 
in more than one dimension. This means that in these cases there is no easy 
way to determine the alternating sums via Morse matchings and the homotopy 
type. In the odd n-case, the bad example shows up for the graph C4.9, where 
-^(64.9) = —2. It is therefore still possible that the independence complexes in 
the other odd n-cases, where Z(C m . n ) is conjectured to be 1, are simple. 

3.4 Graphs on the hexagonal grid 

Take any finite subgraph H' of the hexagonal grid H . As in the case of the square 
grid with even n, there is a natural notion of even and odd vertices given by the 



21 





m 




1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 




3 


(0,2) 


(1,1) 


(1,1) 


(2,2) 


(3,1) 


(3,1) 


(4,2) 


(5,1) 


(5,1) 


(6,2) 


(7,1) 




5 


(1,1) 


(1,1) 


(3,1) 


(3,1) 


(5,1) 


(5,1) 


(7,1) 


(7,1) 


(9,1) 


(9,1) 






7 


(1,1) 


(3,1) 


(5,1) 


(5,1) 


(7,1) 


(9,1) 


(11,1) 












9 


(2,2) 


(3,1) 


(5,1) 


(7,1), (8,3) 


(9,1) 














n 


11 


(3,1) 


(5,1) 


(7,1) 


(9,1) 


















13 


(3,1) 


(5,1) 


(9,1) 




















15 


(4,2) 


(7,1) 


(11,1) 




















17 


(5,1) 


(7,1) 
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(5,1) 


(9,1) 





















Table 2. Homology of I{C m ^ n ) for odd n. 



Z{C 


n,n ) 


m 




1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 




2 


-1 


-1 


1 


1 


-1 


-1 


1 


1 


-1 


-1 


1 




4 


-1 


3 


-3 


5 


-5 


7 


-7 


9 


-9 


11 


-11 




6 


2 


-1 


1 


4 


-1 


-1 


i 


1 


-1 


2 


1 




8 


-1 


3 


5 


5 


3 


7 


1 


1 


-1 


3 


-3 




10 


-1 


-1 


1 


1 


9 


-1 


1 


1 


-11 


-1 


1 


n 


12 


2 


3 


-3 


8 


-5 


7 


8 


9 


-9 


14 


-11 




14 


-1 


-1 


1 


1 


-1 


13 


1 


1 


13 


-1 


15 




16 


-1 


3 


5 


5 


3 


7 


1 


33 


-1 


3 


13 




18 


2 


-1 


1 


4 


-1 


-1 


22 


1 


-1 


38 


1 




20 


-1 


3 


-3 


5 


5 


7 


-7 


9 


41 


11 


-11 




22 


-1 


-1 


1 


1 


-1 


-1 


1 


23 


-1 


-1 


89 




24 


2 


3 


5 


8 


3 


7 


16 


1 


-1 


78 


-3 



Table 3. Z(C m ^ n ) for small m and n, when n is even. Taken from Jonsson [Sj. 
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parity of the vertices. Let G be the graph H' , possibly with some identifications 
made on vertices of the same parity and let O be the set of odd vertices in 
G, where the parity is inherited from H. Then, Aq is a simplex of dimension 
|V(G)| — |0(G)| — 1 and thus contractible. The complex i~o is determined in the 
same way as was done for the square grid with even n. For each wGOwe have a 
maximal face {v g V(G) \0\v& Nbd(u)} € Jo- 
Table |3] shows the homology of the complex I(C^ n ) for small m and n. As 
before, an entry (k,d) means that Hk(I(C^ n )) = Z d and that the homology 
vanishes in all other dimensions. As noted in the square grid case, the entries for 
n = 2 could be derived, for arbitrary m from Proposition ^. 61 We note with interest 
that there seems to be a relation between the dimensions in which homology 
appears for Gf^ and Cj I _ 1 but the data is quite limited. Table [5] shows some 
of the first values of the alternating sum Z{C^ „). They have been generated using 
transfer matrices. We have not found any nice factorisation of these polynomials 
and they seem to have mostly roots which are not roots of unity. 





m 




2 


3 


4 


5 


6 


7 


8 




2 


(2,3) 


(3,5) 


(4,8) 


(5,13) 


(6,21) 


(7,34) 


(8,55) 




3 


(4,5) 


(6,7) 


(7,3) 


(9,22) 


(11,23) 


(12,24) 






4 


(6,3) 


(7,4), (9,1) 


(10,8) 


(11,8), (13,5) 








n 


5 


(7,6) 


(10,11) 














6 


(9,15) 


(11,4), (13,13) 














7 


(11,8) 
















<s 


(12,19) 















Table 4. Homology of I(C^ n ). 



Z(C 


H ) 

Tt,n ) 


m 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 




l 





1 


-1 





1 


-1 





1 


-1 





1 


-1 




2 


2 


-3 


5 


-8 


13 


-21 


34 


-55 


89 


-144 


233 


-377 




3 





-5 


-7 


3 


22 


23 


-24 


-92 


-67 


141 


367 


152 




4 


2 


-3 


5 


-8 


13 


-21 


34 


-55 


89 


-144 


233 


-377 


n 


5 





6 


-11 


-5 


51 


-76 


-60 


416 


-536 


-655 


3351 


-3646 




6 


2 


15 


17 


55 


160 


231 


886 


1664 


3947 


11121 


21065 


59296 




7 





8 


-15 


-35 


57 


34 


-42 


687 


20 


-4207 


-2379 


3611 




8 


2 


-19 


37 


-88 


533 


-725 


3466 


-11927 


21417 


-105552 


273881 


-682665 




9 





-41 


-43 


183 


958 


941 


-9924 


-22943 


19265 


289806 


587437 


-1949599 



Table 5. Z(C^ n ) for small m and n. 



In the case of the hexagonal lattice, we have also calculated the homology for 
the toroidal case, i.e, for the complex I(T r ^ n ). Table [5] shows the homology for 
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small cases. Table [7] lists Z{C^ nn ) for some small values of m and n calculated 
using the transfer matrices. For m < 4, the characteristic polynomials P(m) of 
these matrices factor reasonably well After removing powers of i, the polynomials 
are given by 

P(l) = 1 + t + t 2 
P(2) = 1 + t-t 2 

P(3) = (1 - £ + t 2 )(l - i + 2< 2 + t 3 ) 
P(4) = (1 + t - i 2 )(l - t) 3 (l + i) 3 (l + i 2 ) 

For larger values of m, the polynomials do not factor as nicely, and seem to have 
mostly roots which are not roots of unity. 





m 




2 


3 


4 


5 


6 


7 


8 




2 


(1,3) 


(2,4) 


(3,7) 


(4,11) 


(5,18) 


(6,29) 


(7,47) 


n 


3 


(2,4) 


(4,10) 


(6,4) 


(7,17) 


(9,32) 


(10,1), (11,3) 


(12,76) 




4 


(3,7) 


(6,4) 


(7,15) 


(9,1), (10,12) 


(11,20), (12,1), (13,3) 





Table 6. Homology of I(T^ n ). 





m 


1 


2 


3 


4 


5 


6 


7 


n 


1 

2 
3 
4 
5 
6 
7 
8 


-1 


-1 


2 


-1 


-1 


2 


-1 


-1 


3 


-4 


7 


-11 


18 


-29 


2 


-4 


-10 


-4 


17 


32 


2 


-1 


7 


-4 


15 


-11 


22 


-29 


-1 


-11 


17 


-11 


-51 


127 


-36 


2 


18 


32 


22 


127 


192 


394 


-1 


-29 


2 


-29 


-36 


394 


552 


-1 


47 


-76 


55 


-411 


1478 


83 



Table 7. Z{T% n ) for small m and n. 



4 A second matching 

In this section we introduce an alternative matching on the face poset of the com- 
plex I(C m . n ). While in Section the matching presented was always acyclic, we 
will make no such claims here. This means that we can not rely on the discrete 
Morse theory for our results. Still, we will be able to derive results on the al- 
ternating sum Z(C m ,„). A similar idea was used with great success on the hard 
square model on a torus in [5]. 
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Let M be any matching on the underlying graph of P(J(C TOj „)). Then, for all 
O, t) G M, we have (-1) i<t| + (-l) |r| = and therefore, 

o-e/(C„,„) <?ex 

where we let X denote the set of unmatched faces in I(C m>n ). Thus, we can 
calculate Z(C m ^ n ) using only the faces in the set X. 

In order to describe our matching we need to introduce some notation and 
terminology. Let a be a face in I(C m , n )- Define 

n(a) :={(l,j)€a}. 

Let 7r(er) = {x±, . . . ,Xk} be ordered by the second component of the Xi. For 
1 < i < k, wc define the interval Ni of a at Xi — (1, ji) to be the set 

N t = + 1), + 2), . . . , (1, jt)}, 

where the addition is modulo n in the second component, and the indices and 
k are identified. Let N(a) = {N%, . . . , Nk} be the set of all intervals. We call the 
sequence [|^Vi|, ■ ■ • , \Nk\] the signature of a and wc identify signatures which arc 
identical up to cyclic rotations. We distinguish between even and odd intervals 
by the parity of their sizes and between even and odd positions in an interval 
Ni by saying that the element x% is in an even (odd) position if Ni is an even 
(odd) interval, and extending this notion of parity to the rest of the interval. Let 
7r(cr) = 7r (cr) U 7r e (cr), where n is the set of particles in odd position, and 7r e (cr) 
the set of particles in even position. 




(1,1) (1,4) (1,7) 

Fig. 13. A face a G C 3 ,i 3 . 



Example 4-1- Let a G (^3^3 be the configuration in Figure [13] Here 7r(<r) = 
n (a) = {(1,1), (1,4), (1,7)} and the intervals N(a) = {Ni,N 2 ,N 3 } are given 
by 

m = {(1, 8), (1, 9), (1, 10), (1, 11), (1, 12), (1, 13), (1, 1)}, 
N 2 = {(1,2), (1,3), (1,4)}, and 
tf 3 = {(1,5), (1,6), (1,7)}. 

In the interval A^, the positions (1,5) and (1,7) are odd, while in A^i, the odd 
positions are (1,8), (1,10), (1,12) and (1,1). The signature of a is [7,3,3] which 
is the same as [3, 7, 3] or [3, 3, 7]. 
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We start by defining a closure operator a i— ► a on the poset P(I(C m , n ))- The 
idea is to use the equivalence relation defined by a ~ t if and only if <r = t to 
group faces in equivalence classes. By constructing matchings on the non-singleton 
classes of ~ we can then restrict our study to a particular collection of singleton 
classes. As before, let n(a) = {x\, . . . ,Xk}- The face a is defined as follows. 

— a C a. 

— For each i = 1, . . . , k, if Xi is in an even position and if Xi is free in a, then 
Xi G a. 

Let ~ be the equivalence relation described above and let X a be the equiva- 
lence class of a with respect to ~. That is, 

X a = {re I{C m , n ) | t 

Let a be the face from Example 14. II Then, a = a +(1,9) and X a = {a, a + (1, 9}. 

Lemma 4.1 Let er € I(C m>n ). 

— Let x be a free and even position in an interval N G N(a). Then, a ^ a + x, 
ir (a + x) = 7r (er), and 7r e (<7 + x) = Tr e (a) + x. 

— Assume that \N(a)\ > 1. Letxi G 7r e (er). Then, a ~ cr— a;,, 7r ((T— ajj) = 7r (cr), 
and 7r e (cr — Xi) = 7r e (cr) — x^. 

Proof. The first part follows from the observation that adding the element x to 
a divides the interval TV, into one even interval, with last position in x, and one 
interval of the same parity as Ni, given by the position xi. 

The second part holds since gluing an even interval with its following interval 
(which is necessarily different due to the lower bound condition on N(o~)) creates 
an interval of the same parity as the latter one. The parity of the positions in the 
new interval are the same as in the old intervals and the vacated position Xi is 
free and even in a — Xi, which shows that a ~ a — Xi . □ 

We partition the faces a G I(C m , n ) into three sets P\, P 2 and P 3 , depending 
on the signature and the size of X a . 

— Pi is the set of a for which \X a \ > 1 and for which there is at least one odd 
interval in N(a). 

— Pi is the set of a for which \X a \ > 1 and all the intervals in N(a) are even. 

— P$ is the set of a for which X„ is a singleton. 

Note that Lemma [4.11 implies that if a G Pi, then r G Pi for all r G X a , i G 
{1,2,3}. 

We can use Lemma |4~T1 to describe the sets X a more explicitly. The idea is to 
use the first part to move from a to a. In the process, we add only even intervals. 
Then, we can use the second part to move from a to any r G X a by doing the 
reverse operation of gluing an even interval with its following interval. 

Corollary 4.2 Assume that a € Pi. Then, r G X a if and only if r can be 
obtained from a by removing a subset of ir e (a). 

Corollary 4.3 Assume that a G P%. Then, r G X a if and only if r can be 
obtained from a by removing a proper subset of ir e (a). 
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We can therefore describe the structure of the equivalence classes in the fol- 
lowing way. 

U°\S | S C7r e (a)} if ere Pi, 
X a = U3\S \ S £ tt(c?)} ifaeP 2 , (9) 
{{a} if a G Pa- 

Using the description of X a in ©, we can show that the equivalence classes 
in Pi and P2 behave well. These results are given in Lemma 14.41 and Lemma 14.51 

Lemma 4.4 Let a € Pi. Then, 

£(-l)M=0. 

Proof. Since \X a \ > 1, we know that ir e (a) is non-empty. In particular, we can 
choose an x% € 7r e (<x) and match r — Xi with t + Xi for each r € □ 

Lemma 4.5 Let cr G P2. Then, 

(-i) |t| = (-i)i°'\ 7r ( <7 )i+ 1 . 

Proof. Note that for all cr G P2, we have 7r e (cr) = 7r(er). According to Q, the set 
X CT consists of the faces a \ S where S $1 7r(<7). Therefore, 



(_1)M= Y, (-l)W -(-l)l 5 W ? )l = 0-(-l) 



<t\tt(ct)| 



tEX, \ T=ff\S,SC7r(S 



□ 

Lemma 14.51 needs a bit more work to use than Lemma 14.41 In order to find 
the sum of (— l)' cr ! over all a G P2, we must use Lemma 14.51 on each equivalence 
class in P2. Luckily, there is an easy way to describe all equivalence classes. For 
cr G P2, let Y a = {t G P2 I t \ 7t(t) = cr \ 7r(cr)}. That is, Y a is the set of faces in 
P2 which are equal to a on rows 2 through m. Then, since all faces in P2 contain 
only even intervals, Y„ can be partitioned into two sets, 

y; 1 = {r G Y a I tt(t) C {(1, 1), (1, 3), .... (1, n - 1)}} 

and 

y CT 2 = {r G F CT | tt(t) C {(l, 2), (l, 4), . . . , (1, n)}} 

Now, the set of non-empty Y# and Yj 2 equals the set of equivalence classes in P2. 
This observation leads to the following lemma. 



Lemma 4.6 



£(-l)M=-2 J2 (- 1 )' T| + 

<^6-P2 TGJ(C m _l,„) 



2 .(_l)mn/4 if mis even, 
if m is odd. 
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Proof. Following the previous discussion, we reorder the sum to be over the sets 
Y a and since Y a is uniquely determined by rows 2 through m of ex, we will sum 
over faces [i € C m _i >ra , abusing notation slightly to write = Y a when [i equals 
a \ 7r(cr) as sets. 

EH) H = EEH) M = E ( E (-i) |T| + E 

creP-2 Ya reY„ M6C m _i,„ \reY^ tEY? 

When Y^ and Y^ in the inner sum are non-empty, we can use Lemma 14.51 The 
set Y* is empty precisely when the first row of [i (which is the second row of 
the corresponding a) has particles in positions (1, 1), (1, 3), ... , (1, n — 1). Let X\ 
be the set of such /i and let Xi be the set of /i which has particles in positions 
(1, 2), (1, 4), . . . , (1, n). This gives us the following expression. 

E (-i) H = E 2 • - E (-!) |T|+1 - E (- 1 ) |T|+1 

aeP2 (i€C m _i,„ rgXi t£X 2 

The lemma follows by constructing a matching tree for the set X\ (or equivalently, 
for X2). This set can be perfectly matched when m is odd and can be matched 
except for a single face with mn/4 particles when m is even. □ 

It remains to describe the faces a which are in the set P3. In order to do this, 
we will further divide this set into three disjoint sets, Qi, Q2 and Q3. 

— Qi is the set of faces a with signature [3, 3, . . . , 3]. 

— Q2 is the set of faces a for which all the intervals are of odd length, at least 
one of length greater than 3, and such that all even positions in the intervals 
are non-free. 

— Qs is the set of faces a with signature 0. 

These sets partition the remaining faces since all a with both odd and even 
intervals are in Pi and all a with only even intervals are in Pi- 

Lemma 4.7 

f0 ifm= 3 0, 

E(-1) H = 3-(-l)«/ 3 ifm= 3 l, 

[3 ifm= 3 2. 

Proof. Note that the faces with signature [3, 3, . . . , 3] fall into three equivalence 
classes of ~. Let C' mn = C m; „ \ {(1, 3i + 1) | i = 0, . . . , ra/3 — 1}. We will first 
determine the homotopy type of I(C' mn ). We then obtain the desired sum by 
determining 

E (-i) H+n/3 

and multiply this by 3. 

Figure [14] shows a sketch of the matching tree construction. If we split on v\, 
the subtree where v\ belongs to all faces will be perfectly matched on V2- We can 
therefore assume that v\ is empty, and then match on 113. In the second graph, 
these vertices have been fixed. We now have two vertices of degree 1, which we 
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V2 V3 




Fig. 14. Matching tree for faces a with 7r(<r) = {(1, 3i + 1) | i = 0, . . . , n/3 - 1}. 



can match in any order to obtain the third graph. Note that after removing the 
first three rows and a total of 2n/3 particles, the third graph is identical to the 
graph in the root. This provides a recursive relation. The base cases are when 
m = 1, m = 2 and m = 3. The first two of these leave a single configuration with 
n/3 and 2n/3 particles, respectively. For m = 3, the entire complex is matched. 

!• if m =3 0, 

5 2(rn+l/2W9-l jf m =3 1, 
S2(rn+l)n/9-l if m =3 2. 

The lemma follows. □ 
Conjecture J^.l. Let n be odd. Then, for a e Q2, 

E (-i) |r| = o- 

r :7r(r)— 7r(cr) 

We can prove the following, which shows a special case of Conjecture 14.11 but 
also holds for even n. 

Lemma 4.8 Let a g Q2 and assume that for all intervals iVj G N(a) 7 it holds 
that \Ni\ > 5. TTiera, 

E (-i) M = o. 

r :7r(r )— 7r(cr) 



Proof. Figure [15] shows the fixed vertices in the lower part of a typical config- 
uration in C>2. In this particular example, n — 12 and 7r(cr) = {(1, 5), (1, 12)}. 
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Fig. 15. Configurations from Q2 with n = 12 and 7r(cr) = {5, 12}. 



The same graph, when the fixed vertices are removed, is shown in the root of 
Figure 1161 We will now construct a matching tree on this graph and show that 
we can match all configurations t, with 7t(t) = 7r(er) perfectly. 

First, we match on u\ and V\. By doing this, we obtain two new vertices of 
degree 1, namely ui and i>2. We also have the vertex V3 which had degree 1 already 
in the root node. We continue by matching on U2,i>3 and V2 (in any order), and 
obtain the final graph in Figure [TH Removing all fixed vertices in this graph 
yields, after a rotation of the cylinder, the graph in the root node of height 2 less. 
This is the recursive step. It is easy to see that the base cases, for m = 2 and 
m = 3 are perfectly matched. When m = 1, the set Q2 is empty. 

We argue that this construction is extendable to any 7r(er), where a g Q 2 - Take 
such a a and an interval Ni € N(a). On the second row, there are (|iVj| — 3)/2 
forced particles, due to the fact that all even positions in Ni must be non-free. 
Let yi be the last such particle, which always exists since |iVj| > 5. Then, for each 
i, yi and Xi forces a degree-1 vertex. In the example these degree-1 vertices are u\ 
and v\. We start by matching on these vertices. In the second step, we note that 
for each interval, there will now be (|iVj| — 3)/2 vertices of degree 1 on row 3. In 
the example these vertices are ui for the interval {(1, 1), . . . , (1, 5)} and «2 and ^3 
for the interval {(1, 6), . . . , (1, 12)}. As in the example, when m = 2 and m = 3 
everything can be perfectly matched. Therefore, the number of intervals does not 
matter, as long as they are all odd, and all of size greater than or equal to 5. □ 



The final set Q3 will function as an inductive step, since we have: 

E (- 1 ) kl = E ( 10 ) 

o-eQs ere/(C m _i,„) 



Proposition 4.9 Conjecture \4-l\ implies Coniecture \2. 1\ 



Proof. Let n be odd. We know from Lemma I4T41 that the alternating sum over Pi 
is 0. Since m is odd, we further know that P2 is empty. For P3, we use Lemma FTT1 
for Qi and Conjecture 14.11 for Qi- Finally, as noted in (|10[) . we get an inductive 
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Fig. 16. Matching tree recursion for the configurations from Figure [TH 
step from Q3. 

z(c m , n ) = J2 (-i) w + E (- 1 ) M + E = 

aePi o-eP 2 creP 3 

= 0+0+ E(-!) H + Eh h+ Et- 1 ) 1 ^ 

o-eQi ceQ 2 o-eQ3 

!0 if to =3 and 3|n, 

—3 if to = 3 1 and 3|n, 

3 if to =3 2 and 3|n. 

The base case is Z(Ci. n ) which is the alternating sum of the independence complex 
on an n-cycle. Its values are given by the following expression. 

f-2 ifn= 3 0, 
Z(C hn ) = \ (-l)(«-D/3 ifjlE3lj 
[(_l)(»+l)/3 if n = 3 2. 

Since n is odd, this can be simplified to Z(Ci. n ) = —2 if 3|n and Z(Ci. n ) = 1 
otherwise. Together with the inductive step, this shows the proposition. □ 

When n = 3, the set Q2 is empty and for n = 5, 7 and 9, Lemma l4~8l covers all 
the cases of Conjecture 14. II Therefore, Conjecture 12. II holds for these values of n. 

Conjecture 14. II does not hold for all even n. A counterexample is given in Fig- 
ure [TT1 where m = 3, n = 8 and 7t(ct) = {(1, 3), (1, 8)}. A straightforward applica- 
tion of the matching tree construction produces a matching on the faces t with 
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7t(t) = 7t((t) which leaves a single face with six particles. Thus, X) r -7r(r)={3 8}(~1)' T ' = 
1 in this case. 




Fig. 17. Counter-example to Conjecture |4J] when n is even. 

Still, we can determine Z(C m . n ) for some small, even n. When n = 6, both P\ 
and Q2 are empty. Lemma l4~6l Lemma 14771 and (fTUf gives 

(-l) m /2 m= 6 
m = 6 1 
(-l) m /2 m= 6 2 
m = 6 3 
(-l) m /2 m = 6 4 

m = 6 5 

Therefore, Z(C m fi) is the following periodic sequence, starting from m = 1. The 
period is 12. 

2,-1,1,4,-1,-1,4,1,-1,2,1,1,2,-1,1,... 



E< 


-1)1^1 


m 




1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 




2 


























4 


























6 


























8 






8 




8 




8 


-8 


8 


-8 


8 




10 










10 








-10 






n 


12 














12 












14 












14 






14 




14 




16 






8 




8 




8 


24 


8 


-8 


24 




18 














18 






36 






20 










10 








50 








22 
















22 






88 




24 






8 




8 




20 


-8 


8 


64 


8 



Table 8. Alternating sums over Q2 for small m and n. 



In Tabled] we have taken the values of Z(C m ^ n ) from Table [3] and removed 
the part which can be explained using Lemma l4.6[ Lemma l4.7l and (|10[) similarly 
to what we did in the previous example. What remains in the table is therefore 
the alternating sums over the faces in Q2, or alternatively over those faces in Q2 



E (-d h = 



E 

cre/(C m _i, n ) 



(-l) |CT| + < 



+ 
3 + 2 
3 + 
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which are not matched by Lemma 14.81 These are still unexplained, but several 
interesting patterns seem to emerge. 



5 Transfer matrices 

In [3J, Bousquet-Melou, Linusson and Nevo determined the partition function of 
the hard particle model on a parallelogram shaped region P m ^ n of the square grid. 
We describe it by letting (a, 6) <G V(P m _ n ) if 



Figure [T8l shows the region Ps,6- The skewed border allowed them to find a 
Morse-matching for the complex I(P m ,n) and determine its homotopy type. They 
further realised that this could be used to produce parts of the spectrum of the 
transfer matrix for the cylinder C mi „. In this section, we will review this idea and 
show how we can complete the case for m = 4. We also find the missing eigenvalue 
when m = 6. Unfortunately, we were not able to generalise our findings to new, 
unknown spectra. 

We start by defining the two transfer matrices of interest. For a fixed m, the 
matrix T(m) describes the configurations on the region P mi „ and T'(m) describes 
the configurations on the region S m ^ n . Throughout this section, i will denote the 
square root of —1. 

Let T(m) be the 2 m x 2 m -matrix, indexed by subsets of [m] on both rows and 
columns, and defined as follows. 



where B' = { j + 1 | j <E B n [m - 1]}. 

Let T'(m) be the F m+ \ x _F! m+1 -matrix, indexed by independent sets on a path 
with m vertices, and defined as follows. 



1 < a < m and 1 — a<6<l — a + n. 



(1,1) 



Fig. 18. The graph P 5j6 . 




i\A\+\B\ if AnB = AnB > = 

otherwise, 



(11) 




(12) 
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We here follow [5] and in contrast to the example in Section 12.41 distribute the 
weights of the entries in row A, column B and row B, column A evenly. 

As was illustrated in Section [53 we can determine the values of Z(C m ,n) by 
taking the trace of the nth power of either one of these matrices. Therefore, the 
following relation between the generating functions of Z(C m ,n) holds. 



tr (1 - tTim))- 1 -2 m = J2 Z(P m ,n)t n = tr (1 - tT'(m)) * - F m+1 (13) 

n>l 



We conclude that the spectra of T(m) and T'(m) are equal, apart from the mul- 
tiplicity of the null eigenvalue due to the difference in size. 

Let A,BC [m] and assume that we prescribe the particles given by A and B 
on the left and right boundary of P m ^ n . That is, when j S A, then the jth particle 
from the top on the left boundary is present, and otherwise it is missing in each 
configuration and similarly for B on the right boundary. This induces a subgraph 
of P m ,n which we will denote by P rnj „(A, B). 

Let GA,B(t) denote the element is row A, column B of (1 — tT(m)) -1 . Then, 
GA,B(t) is the generating function for configurations on P m n (A, B), weighted as 
given by (fTTj) and has the following expression. 

G A , B {t) = 6 a ,b + H) |A|+|S| Z(Pm,n+i(A, B))t n (14) 

n>l 

We can then split the left-hand side of (fT3"| into 

tr(l -tT(m))- 1 = J2 G A,A(t)- (15) 

AC[m] 

The observation made in [3] was that parts of the spectra can be recovered 
by finding GA,A{t) for some A and that this can be done in the model for P m ^ n 
which is more easily determined than that of C m>n . 

The characteristic polynomial of T'(4) is determined both in [3] and [H] and 
is given by (1 - t + t 2 )(l - t 2 )(l - i 4 ). Let G A N) B (t) denote the first N terms in 
the expansion of (fT4|) . By finding a Morse-matching for the region in Figure IT8| 
Bousquct-Mclou, Linusson and Nevo effectively calculated G$$(t). When m = 4, 
they found that 

G$m = 1 _] +t 2 - ( 16 ) 

This determined 

This recovers one of the factors of T'(4). Using this case, and a recursion for 
Z (P mjn ({2, 3}, {2, 3})), they could then determine the function G{2.3},{2,3}(i)- 

G {a ,3 } , { 2,3}(*) = (i+^^t*) W 

We will now determine the remaining cases for m = 4. We start with the 
following lemma. 

Lemma 5.1 Assume that j — l,j,j + 1 G A for some j and let B = A \ {j}. 

Then, 

G A>A (t)+G B ,B( t ) = 2 - 
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Fig. 19. Two regions differing only in the position j. 



Proof. Figure [H] illustrates that the regions remaining after fixing the boundaries 
to A and B respectively are equal. Since the addition of two elements does not 
change the parity, we have Z (P m ^ n (A, A)) = Z (P m ,n{B, B)). Since furthermore 
(_ _ — ( — all terms except the constants cancel and we have 

G A! A(t)+G B! B(t)=dA,A + 6 B ,B = 2. □ 

Using Lemma 15.11 we can avoid calculating certain sets by choosing a partial 
matching on 2^ which satisfies the condition in the lemma. We will in this way 
match {1, 2, 3} with {1, 3}, {2, 3, 4} with {2, 4} and {1, 2, 3, 4} with {1, 3, 4}. This 
leaves only {1,2,4} unmatched among the potential candidates for Lemma f5.1l 
We can produce its generating function by noting that G[^^(t) = 1 + t 2 G®$(t). 
Lemma l5~T1 applied to A = [4] and i = 3 gives us 

1 - t 

T+t 

The sets that remain are the singleton sets, {1,2}, {3,4} and {1,4}. By sym- 
metry, we can restrict our calculations to the sets {3}, {4}, {1,2} and {1,4}. We 
present the calculations for {4} and {1,4}. 

We start with the case A = and B = (3, where j3 is an arbitrary subset of 
[to] . When j > 4, we have the recursive relation Z (P 4 j (0, /?)) = —Z (P 4 j-3 (0, /?)). 

T ®,/3 



G{1,2,4},{1,2,4} W = 1 - * G$$(t) - l t+t 2 - ( 18 ) 



Thus, Gtft t p(t) — G^Jt)— i 3 G0 j( g(i) which gives us the following general expression. 



- YT^ (19) 



The recursive relation which was used to produce (|19j) . and those that follow, 
are obtained using simple matching trees. We present some, but not all trees here. 
Figure [2171 illustrates a matching tree when A = {4}. This gives us the following 
recursive relation for the case A = B = {4}. 

Z (P 4J ({4}, {4})) = Z (P 4 ,,- 4 ({4}, {4})) + Z (P 4j - 4 (0, {4})) (20) 

For the second term, we calculate G 4 | 4 |(t) = —it, which, using (|19|) . gives 
G®.{4}(t) = — it/ (1 + t 3 ). Plugging this term into (j2TI|) gives the expression 

G {4} , {4} (t) = G\% {4} (t) +t 4 (G {4}i{4} (f) - 1) +it 4 G $M} {t), (21) 

which, after calculating { 4 }(^) = 1 + ^ 4 : gives the final function for A = B = 
{4}, and by symmetry for A = B = {1}. 

1 + t 3 + t 5 

G{4}.{4}(t) = G { 1 } , {!}(<) = +f 3)(! _ ^4) ( 22 ) 
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Fig. 20. The large circles indicate match nodes and the large triangle indicate 
split nodes. The numbers give the order in which the nodes are matched or split. 
At 5, the two branches of the recursion are obtained. 



For A = B = {1, 4}, we have the recursion 

Z(P^(A,A)) = Z(P 4 j-i(0,A)) + Z(P 4 j-4(A,A). (23) 
We proceed by summing over j and calculating the base cases. 

G{l,A},{l,4.}(t) 

= Gf^MM} W - t (G(»,{M}(t) ~ G 4 {M} (*)) + (%,4},{i,4}W - 1) 

= (W+t 4 + T ^-t 2 -^/(l-t 4 ) (24) 

In the last equality, we have calculated the polynomials 

G {?,4},{i,4}(*) = 1 + <2 + * 4 and G um® = 

and used the expression for 00/1,4} (i) in (fT9|) . It turns out that Gs\^\ ,{1,4} (t) = 

G{2,3},{2,3}(*)- 

The final two cases, when A = B = {3} and A = B = {1,2}, are determined 
similarly, and we present them without any explicit calculations. 

1 + t 3 + t 5 

G{2},{2}(*) = G {3},{3}W = (! + f 3)(! _ t 4) ( 25 ) 

Note that (|25|) implies that all the singleton sets have the same generating func- 
tion. 

t b 

G{l,2},{l,2}(i) = G{3,4},{3,4}(*) = 1 - ( 1 + t 3)(! _ t 4) - ( 26 ) 

We note that while all factors of the characteristic polynomial show up in 
the individual terms, we can not deduce their multiplicities until we have taken 
the sum over all sets. Expanding this sum into partial fractions, we arrive at the 
following expression. 



GA,A(t) = G j(t) + G{ 1i2 ,4},{1,2,4}(0 + 6+ 

AC [4] 

+ 4G{4},{4}(£) + 2G{ 2 ,3},{2,3}(*) + 2G{x,2},{l,2}{t) 

_ 1 1 2 2 1 1 

+ l-it + l + it + l + t + l-t + l-at + l-at 1 
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where a = (1 + i\/3)/2 and a are the two roots of 1 — t + t 2 . As expected, we 
recover the roots and multiplicities of the characteristic polynomial of T(4), which 
apart from the 8 zeroes equals the characteristic polynomial of T"(4). 

Next, we determine the factor of the characteristic polynomial for m — 6 which 
does not follow from the case A = B = 0. The characteristic polynomial for m = 6 
is given by (1 - i 4 ) 2 (l - t u )/(l + t) [5|5] and the function G , (i) follows from 
the general result in [3] . 

i / ^ -t 12 \ 

This leaves the factor 1 — t 4 unexplained. We show that this factor appears in the 
generating function Ga,a(*) for A = {2,3,6}. 

As in the previous case, we start by deriving an expression for G$ p(t) for an 
arbitrary (3. The matching leading up to the recursion contains more steps, but 
is straightforward using matching trees. The final formula is 

G *dt) = y^t- (28) 

We will also make use of the recursive relation 

Z (P 6J ({2}, {2, 3, 6})) = -Z (P 6>i -6(0, {2, 3, 6})) , (29) 
from which we derive that 

(7) 6 
G{2},{2,3,6} W = G {2} {2 3 6} (*) + G 0,{2,3,6} (*) 

= iW (zffiJ^!i) _ tOg+fl . (so) 

In the last equality we have used that G^ 1 ^ 3 6 y(t) = —i(t + t i + t 12 ). The 
main recursion for A = {2, 3, 6} is given by the following relation. 

Z (P 6 j(A, A)) = Z (P 6>j -i(A, A)) - Z (P 6J _ 1 (0, A)) - Z (P 6j _ 4 ({2}, A)) (31) 

The matching tree needed to produce this relation is shown in Figure 1211 The 
terms Z (P 6)J -_i(0, A)) and Z (P 6iJ _ 4 ({2}, A)) in the right-hand side of flU} are 
obtained in nodes 1 and 4, respectively. After summation over i we have 

G{2,3,6},{2,3,6}M = G {2,S,G} ,{2,3,6}^) + ^ (^{2, 3, 6}, {2, 3,6} (*) _ l) 

- it (G 0:{2 ,3,6} (*) ~ G 4 {2,3,6} (*)) + ^{2}, {2, 3, 6} (*) (32) 

Calculations of initial terms in the relevant series provides the polynomials 

^{2,3, 6}, {2, 3,6} (*) = 1 — * + t an d G j 2 >3i6 }(*) = ~ **) 

so (f3"2"|) evaluates to 



r f.s (l+t-^-^-f+t^+t' 2 )(l-t) 
^{2,3,6}, {2,3,6}W - (1 - < 14 )(1 - i 4 ) ' 
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Fig. 21. A matching for P 6j ({2, 3, 6}, {2, 3, 6}). 



This recovers the missing factor 1 — t A . 
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